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We study the superconducting transition temperature, T c , of an array of fluctuating, maximally 
gapped interacting ladders, embedded in a disordered plane. Renormalization group analysis indi- 
cates that geometrical fluctuations mitigate the suppression of the intra-ladder pairing correlations 
by the disorder. In addition, the fluctuations enhance the Josephson tunneling between the ladders. 
Both effects lead to an increase of T c in the Josephson coupled array. These findings may be relevant 
to the understanding of the 1/8 anomaly in the lanthanum based high-temperature superconductors. 
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I. INTRODUCTION 

Both perturbative renormalization group (RG) 
studies^ and strong coupling numerical calculations 2 
have demonstrated that doped t — J and Hubbard 
ladders typically exhibit a maximally gapped phase, 
with a single remaining gapless charge mode protected 
by symmetry 3 . Such a phase possesses substantial 
superconducting pairing correlations with approximate 
rf-wave symmetry. By coupling many ladders together 
via Josephson tunneling, the phases of the different 
Cooper pairs may establish coherence, in which case 
global superconductivity ensues. It has been suggested 
that the underlying mechanism of superconductivity in 
the cuprates is essentially the one described here for the 
quasi-one-dimensional system^. 

Previous studies have shown that Josephson coupled 
arrays of spin-gapped chains and two-leg ladders are par- 
ticularly prone to the deleterious effects of disorder on 
T^i On the other hand, we have demonstrated that ge- 
ometrical fluctuations of a spin-gapless one-dimensional 
system tend to reduce the effects of disorder—. It is there- 
fore interesting to try and marry these two observations, 
and study the interplay between disorder and fluctua- 
tions in a quasi-one-dimensional superconductor. Beyond 
its purely theoretical appeal this question may also be 
relevant to a distinctive feature in the phenomenology of 
the lanthanum based cuprates. 

Among the high-temperature superconductors, the 
lanthanum compounds stand out as the ones with the 
most extensive experimental evidence for the existence of 
quasi-one-dimensional charge and spin stripe order, both 
in their " normal" and superconducting states^^. Exper- 
iments have also found clear correlation between the pres- 
ence of robust static stripe order and the suppression of 
superconductivity in these materials^. For example, in 
Lai.e_ a; Ndo.4Sr x Cu04 the ordering temperature, T m , of 
the incommensurate magnetic structure reaches a max- 
imum around x = 1/8, while the superconducting tem- 
perature, T c , exhibits a substantial dip there. This 1/8 
anomaly was first identified in La2- 2; Ba 2 ;Cu04— , where 
it is especially pronounced. In view of the theoretical 



insights outlined above, it is tempting to associate the 
suppression of T c with enhanced disorder effects due to 
the slowing down of the stripe fluctuations. 

Motivated by such questions we study in this paper 
a model of a fluctuating disordered array of maximally 
gapped n-leg ladders. We find that the fluctuations re- 
duce the disorder induced degradation of the pairing cor- 
relations on each one-dimensional system, although they 
typically do not change the critical value of the Luttinger 
interaction parameter, K* , above which infinitesimal dis- 
order is irrelevant. The latter is a decreasing function 
of the width n. Moreover, the fluctuations enhance the 
Josephson tunneling between the ladders. Both of these 
effects lead to a higher T c for the array. 



II. THE SINGLE LADDER MODEL 

We begin by considering a single fluctuating n-leg 
ladder. We assume that the fluctuations occur on 
length scales which are long compared to the electronic 
wavelength and thus neglect backscattering from them. 
Specifically, we treat the limit where the ladder oscillates 
rigidly inside a parabolic potential. Its state is given by 
the deviation u(t) from the bottom of the well (taken in 
the following as the a:-axis), and its action in imaginary 
time is 



S u = dr 



-u 



(1) 



where M is the ladder's effective mass. The fluctuations 
are characterized by the oscillation frequency u>q, and by 
their typical amplitude A = I/^/AIluq. 

The non-interacting spectrum of the ladder consists 
of n energy bands. We assume that the filling is such 
that all of them cross the Fermi level. If this is not the 
case, then, in the following, n should be taken as the 
number of bands that do cross Ep . Each of those do 
so at two Fermi points ±kp\ one containing left-moving 
(r = —1), and the other right-moving (r = 1) particles. 
The annihilation operator of an electron in band b, with 
chirality r and spin projection a = ±1 can be expressed 
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using bosonization language as 



ipb,r,a(x) 



v'2 



(2) 

where the Klein factors F(, iT . a ensure the correct anti- 
commutation between the different particle species, and 
a is a short distance cutoff of the order of hp 1 . Phys- 
ically, (frc^ and (j)^ are, respectively, the phases of the 
2kp^ charge and spin-density wave fluctuations in band 

b, while 6 c b ^ is the superconducting phase of the band. 

When the electronic interactions are turned on, one 
finds for generic densities away from half filling that the 
system ends up in a phase where all the spin fields <pi , 
and all the relative charge modes 8c — 8c, acquire 
a gapii£. Consequently, the pronounced fluctuations of 
their conjugated fields di b \ and (f> c — 4> c h \ cause their 
correlations to decay exponentially to zero. In this so 
called C1S0 phase the only gapless excitations are cre- 
ated by the total charge fields 
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Note that if the interaction strength is larger than the 
band-structure gaps, the number of occupied bands may 
change relative to the occupancy of the bare spectrum, 
and n should be adjusted correspondingly. 

The plane containing the ladder is taken to include a 
weak random potential which couples to the electrons, 
and through them to u. This choice is motivated by the 
situation in the cuprates where the "ladders" are defined 
by the sites occupied by holes, and as such interact with 
the disorder only indirectly via the electronic degrees of 
freedom. Consequently, the forward, V n (r), and back- 
ward scattering, V*(r), disorder components couple as 

~ n 

H dis = jd X [K,b( x >yj) + h ib( x >y^] +B -- c -' ( 4 ) 

a,b,j=l 

where yj = jd + u, with d the inter-leg distance, and 
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The exponents in Eq. (|5|) are given by 

Xn ,d^b,r,a) = T - (0<«O - <#»)) 

+ r<r(0to t - <r(0W - (6) 

They vanish for the intra-band forward scattering pro- 
cesses and instead these terms take the form 

KM ~ K!a(r) (n a + —d x ^\ , (7) 



where n a is the average electronic density in band a. 
Assuming a short-range Gaussian disorder, one finds 
that its components obey V^ Ar (x, yj)[V% ^ y (x, yf)]* = 

DB al [&r,r'5a,ai5b,b> + 5 r - r >5 a ,b>5 bta i]5(x - x')S(yj - yf), 

mdVl b!r (x,yj)[V$ !b , y (x,yf)}* = DB^5 ry [5 a , a ,6bM+ 
8a,b'8 b ,a']fi(x — x')6(yj — yf)- Here, D is the disorder 
strength and B^ b are numbers of order unity which de- 
pend on the details of the band structure. 

In the following we focus on the physics at energies of 
the order of the inter-ladder couplings, which we assume 
lie well below the scale A, set by the gaps of the massive 
modes. Under such conditions the latter can be replaced 
by their non-vanishing expectation values, thus yielding 



X ,, 5 (a, b, r, a) = T ^ ] ) ~ ^ - 0^) 

+ c n .(:(a,b,r,cr), (8) 



where c^f (a, b, r, cr) are constants. All the Xv,C contain 
fields conjugated to the massive modes, whose correla- 
tion functions decay exponentially. As a result, except 
for the intra-band forward scattering part (whose effects 
we discuss in the next section) the bare disorder potential 
coupling, Eq. is irrelevant in the RG sense. However, 
higher order terms, generated in the course of the RG by 
integrating over the massive modes, may be relevant£ii£. 
To obtain the structure of these terms, which are the 
ones responsible for the degradation of the pairing cor- 
relations on the ladder, we need to consider all possible 
linear combinations of the Xn,i with integer coefficients. 
The most relevant coupling is that which contains the 
gapless field C + with the smallest coefficient in front of 
it, and no gapped fields. Eq. © indicates that this cou- 
pling is genericallyi^ generated by Ylb=i Xf (&j ^ r i a )i an d 
takes the form 2J2b=i 0° = %\ / n ( t > c+- Consequently, the 
scattering by the effective disorder is described by 



ff-dis 



dx 



V c ii(x,u) 
27ra 



H.c. 



(9) 



where a is of the order of the correlation length v/A of 
the gapped fields, with v being their typical velocity^. 
The leading contrib ution to the effective disorder satisfies 
V eS {x,y)V* s {x>,y>) = D cS S(x - x')5(y - y'), with 



D cS ~ Dp/aAA 2 )"- 1 (A/o;o)("- 1 )/ 2 (a/a) 



n+l 



(10) 



Hence, after averaging over the disorder and the ladder 
dynamics, as described by Eq. ([T]), we are led to consider 
the following action for the CI SO phase of the n-leg ladder 
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S 



dxdr 



-id x 6 c+ d T <p c+ 



V c +K c+ 2 V c+ 2 

^ (aA+) +^—(d x (l)c+) 



c+ 



(27ra) 2 

where the fluctuations kernel is given by^ 

F(t-t') = {S[u(t)-u(t')]) 



(t - t') cos I [(/> c+ (a;, r) - </> c+ (.t, t')] | , 



r 



i 



/2ttX \J\ - e -w \r-r'\ 



(11) 



(12) 



III. RG ANALYSIS 

The RG analysis of this action is most naturally formu- 
lated in terms of the following dimcnsionless parameters 



V = 



nD cS a 



au>o 

V c + ' 



K = K, 



c+ 



V ( 4 , 
4 U 2 



ds- 



(13) 



Note that the parameter K is no longer a pure measure 
of the electronic interactions, but is rather an admixture 
of the interactions and the disorder—. However, it is 
K which controls the asymptotic decay of the various 
correlation functions—. The resulting RG flow equations, 
with respect to the running cutoff a(£) = ae are 



dK 
~~d£ 

dm 

~dl 

d\ 
di ~~ 



1 



V 



2 VT^IP- 

ACT) (coshn? 



dV 
~d£ 



(3 - nK)V, 



1) 



n (l-e- ro ) 3 / 2 ' 
2CD (1 — coshtu + vasmhru) 



(14) 



7)3/2 



with the dimensionless ladder length C = L/2ira scal- 
ing as dC/dl = —C. The details of their derivation are 
presented in Ref. H. 

The RG equations (fl"4"l) imply that infinitesimal disor- 
der is irrelevant for K > K* = 3/n, independently of the 
fluctuations. More important, however, are the effects of 
fluctuations in the presence of finite disorder. First, since 



V oc A" 



-(n-i)/2 



, see Eqs. (fTUf and p3|) . large ampli 



tude and frequency oscillations tend to decrease the ini- 
tial strength of the effective disorder. Secondly, Eq. (fl4|) 
shows that the disorder is more efficient in renormalizing 
K for smaller oscillation frequencies, and that it drives va 
towards this regime. Due to both mechanisms, stronger 
and faster fluctuations cause K to flow more slowly to 
smaller values. Consequently, the ladder's superconduct- 
ing pairing correlations, which decay as x~ 1 ^ nK , are less 
degraded by the disorder. As demonstrated below, these 
stronger pairing correlations translate into an increase of 
T c for the Josephson-coupled array. 



Few comments are appropriate at this point. First, 
note that formally the static ladder limit is reached in 
our model by taking A — > 0, loq — > oo, in order to keep 
the ladder in its ground state inside the well. Secondly, as 
mentioned in the previous section, the intra-band forward 
scattering terms, Eq. ([7]), survive in the gapped phase. 
They are relevant in the RG sense: the dimensionless 
strengths of the parts which couple to the average band- 
density and to its long-wavelength fluctuations scale as 
dT>f y i/d£ = 32?/,! and dDj^jdl = I?/,2, respectively. 
These processes can be shown not to affect the pairing 
correlations in the static limits, or more generally when 
luq — ► oo. However, for finite loq they induce a retarded 
attraction between the electrons via exchange of ladder's 
oscillation quanta. This attraction tends to renormalize 
K towards higher values and enhance pairing, thus re- 
inforcing the effects outlined above. The intra-band for- 
ward scattering processes do contribute to the localiza- 
tion of the ladder since they provide additional channels 
for its coupling to the disorder [beyond the one contained 
in Eq. (fTTj) ]. Consequently, in the last two flow equations 
of set (HU) one should substitute V -> V + V fA + V sa - 
Due to its large scaling dimension X>/,i is particularly 
effective in localizing the ladder. 

Finally, one should also take into account the renor- 
malization of K at energy scales which span the range 
between Ep and A, i.e., during the generation of the gaps 
which define the C1S0 phase. Perturbative RG analysis 
of the single chain problem^' 16 indicates that K renor- 
malizes by an amount of the order of D/vpXA during 
this stage. In the following we assume that D <C ^fAA, 
and ignore this "high energy" renormalization. More- 
over, such a condition is necessary in order to ensure 
that disorder does not destroy the gaps. Since the size of 
the gaps is known to decrease rapidly with the number 
of legs^^iii, the region of applicability of our treatment 
shrinks correspondingly with n. 



IV. COUPLED LADDERS 

Next, we consider an array of coupled CI SO n-leg lad- 
ders. The arguments leading to Eq. ^ imply that the 
most relevant inter-ladder charge-density wave (CDW) 
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coupling, V, is the 2nkp CDW coupling. Its scaling 
dimension, 2 — nK, makes it more relevant than the 
Josephson tunneling only once K < 1/n. Moreover, h2 a , 
which we neglected in this work, is known to suppress V 
exponentially 16 . Inter-ladder q ~ couplings also pro- 
mote J over V^. Owing to these reasons and our inter- 
est in applications to the superconducting cuprates, we 
assume that the Josephson tunneling between adjacent 
ladders dominates over their CDW coupling. Hence we 
analyze the effect of 



Hj 



d.r 



E 



Jab{Ui 



H.c, (15) 



a,b=l 



with A* = ^* i+a ^,-,t + ^i,-4^,+,t the singlet su- 
perconducting order parameter in band a of the i-th 
ladder. The Josephson tunneling amplitude depends, 
exponentially, on the inter-ladder separation, J a b ~ 
Joe~\ s+Ui ~ Uj l^ 7 , where s is the mean spacing of the array 
and 7 depends on the environment between the ladders^. 
Assuming A < s, in order to remain in the quasi-one- 
dimcnsional limit, we obtain for Hj in the CI SO phase 



Hj 



J e («»-«*+l)/7 



■ cos 



J c+) 



_ (16) 

where 1 ! J = 8££ 6=1 Ja 6 (0)(cos(V24 a) )<cos(V24 b) )- 
Its contribution to the RG flow, in terms of J — 
J /2iTy/nv c + 1 and derived assuming 2 ^ J, A/7 -C 1, is 



dK 
~dt 
dm 



= J 2 



r(l 



dl ~ ( 2 tik) J > 



= w- 2mr 2 £J 2 — [1 + I (zu)] , 



7 



(17) 



r 



w 



where Iq^x) are modified Bessel functions. 

By integrating Eqs. (fT4"]) and (fT7|) one can calcu- 
late the scale, £*, at which J grows to be of order 
1, and obtain an estimate for the transition tempera- 
ture of the array T c ~ Ae~ e . Since fluctuations tend 
to increase K, via its renormalization by both T> and 
J : and since K controls the scaling of J, T c likewise 
rises. Another estimate for T c can be derived using 
the inter-ladder mean field theor y 6 ! 7 ' 21 . By replacing 
j e (u t -u ] )/ 1 J = 2J(e( u '-^)/f) = 2Je A2 / 2 T 2 , and 
cos[ v /27^(^ + - 9i + )] -► (cos(y27^6» c+ ))cos( v /27^6' c+ ) 
in Hj, the problem turns into that of a single ladder 
described by Eq. (TTTj) and coupled to an effective field 



Hj 



— / dx- 



W 



2ita)'- 



■ cos i 



(18) 



induced by the neighboring ladders and determined self- 
consistently from W = J (cos ( \J1j n9 c + ) ) . At T c , when 
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FIG. 1: T c of an array of 2-leg ladders as obtained by inte- 
grating the RG equations (| 14f) and (|17|) (solid line), and from 
the mean field theory, Eq. (|20p . (broken line). The results 
are for a system with bare K = 0.7, J = 10" 2 , V = 10" 4 , 
A/7 = 10~ 3 , and £ = 10. 



W is vanishingly small, this condition reads 
1 1 r fi/r. 



J (27ra) 5 



drP(x, r), 



(19) 



with P(a;,r) = (cos[v^7^c+(a:, t)] cos[ v /27^6» c+ (0, 0)]) 
calculated using the action (fTTj) . Approximating the ef- 
fect of temperature by a cutoff length Z<r = w c + /T beyond 
which all correlators decay exponentially to zero, and us- 
ing the scaling properties of P , Eq. (fTTJ|) becomes 6,7 



1 

1 



RdR 



■ exp 



\n(R/a) 



df 



nK(£) 



(20) 



where K(£) is calculated from Eq. |T 

Results for the ratio of T c to the transition temper- 
ature of the pure system, in the case n = 2, are pre- 
sented in Fig. Q] The effect of the fluctuations on T c is 
substantial 22 and is particularly important whenever D 
is more relevant than J, i.e. for K < (l+y/b)/2n. Fig. CD 
was generated for fixed bare T>, ignoring the dependence 
of the effective disorder on wo , as given by Eq. (fT0|) . Tak- 
ing this dependence into account will only make the miti- 
gating consequences of the fluctuations more pronounced. 
Additional amplification of the physics discussed here oc- 
curs when shape deformations of the ladders are taken 
into account^. 



V. DISCUSSION 

Our study was motivated by the physics of the cuprate 
high-temperature superconductors, especially the 214 
compounds, which exhibit both static and fluctuating 
charge and spin stripe orders over a substantial portion 
of their phase diagram^ 5 -. In particular, we were inter- 
ested to gain insight into the anomalous suppression of 
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T c at x — 1/8, which is manifested to various degrees in 
these systems. One should be cautious, however, when 
trying to apply the model studied here to the physical 
materials. To begin with, we have considered the quasi- 
one-dimensional limit of weakly coupled ladders. This 
limit allows for a controlled theoretical treatment of the 
problem, but is an extreme caricature of the actual charge 
modulation observed in the cuprates, whose amplitude is 
difficult to ascertain but which is most likely weak. We 
then assumed smooth long-wavelength fluctuations and 
ignored lattice commensurability effects, which are be- 
lieved to be important in the slowing down of the stripe 
dynamics at x = 1/8. In a related approximation we 
did not include fluctuations of the stripe order towards 
a similar order oriented along the perpendicular direc- 
tion. Recent scanning tunneling spectroscopy^ demon- 
strated domains of perpendicular stripe orientations in 
Bi2Sr2CaCu20s+5. It is conceivable, however, that in 
the low temperature tetragonal phase of La2^ x Ba x Cu04, 
where the in-plane rotational symmetry is broken and 
the x = 1/8 anomaly is observed, such an approxima- 
tion is not too severe. Finally, the striped model pos- 
sesses a large spin-gap, which implies the absence of gap- 
less nodal quasiparticles in the superconducting phase 
(although it exhibits a " d- wave-like" order parameter). 
The same spin-gap means that the model does not con- 
tain any of the low-energy incommensurate spin excita- 
tions, which are the primary experimental evidence for 



stripes. These excitations originate from the intervening 
lightly doped regions between the charged stripes. How- 
ever, apart from acting as effective barriers for tunneling 
between the stripes [as characterized by the phenomeno- 
logical length 7 in Eq. (fT6f ] they have little effect on 
the charge response of the system governing its super- 
conducting susceptibility, which is the focal point of the 
present study. 

Notwithstanding, our model constitutes a controlled 
theoretical laboratory for the study of the interplay be- 
tween the stripe dynamics and the disorder. Provided 
that w is empirically associated with the distance from 
x = 1/8 in the phase diagram of the lanthanum based 
cuprates, then Fig. [1] offers a mechanism for the T c 
anomaly at x = 1/8. It suggests that the above men- 
tioned interplay plays an important role in alleviating, 
at least partially, the disorder-induced suppression of su- 
perconductivity in these compounds. 



Acknowledgments 

It is a pleasure to thank E. Arrigoni and T. Gia- 
marchi for useful discussions. This work was supported 
by the United States - Israel Binational Science Foun- 
dation (grant No. 2004162) and by the Israel Science 
Foundation (grant No. 538/08). 



1 L. Balents and M. P. A. Fisher, Phys. Rev. B 53, 12133 

(1996) ; H. J. Schulz, ibid. 53, R2959 (1996); H.-H. Lin, 
L. Balents and M. P. A. Fisher, ibid. 56, 6569 (1997); 
S. Baruch and D. Orgad, ibid. 71, 184503 (2005). 

2 R. M. Noack, S. R. White, and D. J. Scalapino, Phys. 
Rev. Lett. 73, 882 (1994); C. A. Hayward, D. Poilblanc, 
R. M. Noack, D. J. Scalapino, and W. Hanke, ibid. 75, 926 
(1995); S. R. White and D. J. Scalapino, Phys. Rev. B 57, 
3031 (1998); S. R. White and D. J. Scalapino, ibid. 55, 
R14701 (1997). See also Refs.[lE 

3 V. J. Emery, S. A. Kivelson, and O. Zachar, Phys. Rev. B 
59, 15641 (1999). 

4 E. W. Carlson, V. J. Emery, S. A. Kivelson, and D. Orgad, 
in "Superconductivity: Novel Superconductors", Vol 2, p. 
1225 , edited by K. H. Bennemann and J. B. Ketterson 
(Springer- Verlag 2008). 

5 S. A. Kivelson, I. P. Bindloss, E. Fradkin, V. Oganesyan, 
J. M. Tranquada, A. Kapitulnik, and C. Howald, Rev. 
Mod. Phys. 75, 1201 (2003). 

6 Y. Suzumura and T. Giamarchi, J. Phys. Soc. Jpn. 58, 
1748 (1989). 

7 E. Orignac and T. Giamarchi, Phys. Rev. B 56, 7167 

(1997) . 

8 U. London, T. Giamarchi, and D. Orgad, Phys. Rev. B 73, 
134201 (2006). 

9 J. M. Tranquada, B. J. Sternlieb, J. D. Axe, Y. Naka- 
mura, S. Uchida, Nature (London) 375, 561 (1995); M. Fu- 
jita, H. Goka, K. Yamada, and M. Matsuda, Phys. Rev. 
Lett. 88, 167008 (2002); M. Fujita, H. Goka, K. Yamada, 



J. M. Tranquada, L. P. Regnault, Phys. Rev. B 70, 104517 
(2004). 

J. M. Tranquada, J. D. Axe, N. Ichikawa, A. R. Mood- 
enbaugh, Y. Nakamura, and S. Uchida, Phys. Rev. Lett. 

78, 338 (1997); N. Ichikawa, S. Uchida, J. M. Tranquada, 
T. Niemoller, P. M. Gehring, S.-H. Lee, and J. R. Schnei- 
der, ibid. 85, 1738 (2000) 

A. R. Moodenbaugh, Y. Xu, M. Suenaga, T. J. Folkerts, 

and R. N. Shelton, Phys. Rev. B 38, 4596 (1988). 

E. Arrigoni, B. Brendel, and W. Hanke, Phys. Rev. Lett. 

79, 2297 (1997); E. Arrigoni and S. A. Kivelson, Phys. 
Rev. B 68, 180503(R) (2003). 

In general it can also be constructed using other combi- 
nations. For example, in the 2-leg ladder it follows from 
Xe(l,2,r, ( j) +x«(l,2,r,-a). 

The prefactors in Eqs. ((9} and (|16[) contain combinations 
of products of Klein factors. We assume that we work in 
the subspace of states where their eigenvalue is 1. 
Consequently K and K c + differ in their RG flow. The latter 
actually converges towards a seperatrix at K c+ = 2/n. 
This can be seen by transforming to <f> = \Jn/24> c + and 9 = 
\/2/n(9 c +, in terms of which the action describes spinless 
fermions moving on a single disordered chain, with K = 
(n/2)K c+ . For K c+ — 2/n the fermions are non-interacting 
and the disorder can not change this fact. See Refs. l8lflq 
and I. V. Gornyi, A. D. Mirlin, and D. G. Polyakov, Phys. 
Rev. B 75, 085421 (2007). 

T. Giamarchi and H. J. Schulz, Phys. Rev. B 37, 325 



6 



(1988). 

17 S. Chakravarty, Phys. Rev. Lett. 77, 4446 (1996). 

18 V. J. Emery, E. Fradkin, S. A. Kivelson, and T. C. Luben- 
sky, Phys. Rev. Lett. 85, 2160 (2000); A. Vishwanath and 
D. Carpentier, ibid. 86, 676 (2001). 

19 S. A. Kivelson, E. Fradkin, and V. J. Emery, Nature (Lon- 
don) 393, 550 (1998). 

20 The latter condition allows us to expand the exponential 
dependence of J on m — u i+ i, see Eq. (fi"5jl . 

21 D. J. Scalapino, Y. Imry, and P. Pincus, Phys. Rev. B 11, 



2042 (1975); E. W. Carlson, D. Orgad, S. A. Kivelson, and 

V. J. Emery, ibid. 62, 3422 (2000); E. Arrigoni, E. Fradkin, 

and S. A. Kivelson, ibid. 69, 214519 (2004). 

The fact that T c ~ T CQ for large zu in Fig. 1 is a result of 

the small bare T> f J ratio used in the calculation. 

Y. Kohsaka, C. Taylor, K. Fujita, A. Schmidt, C. Lupien, 

T. Hanaguri, M. Azuma, M. Takano, H. Eisaki, H. Takagi, 

S. Uchida, and J. C. Davis, Nature (London) 315, 1380 

(2007). 



